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Abstract 

We discuss the consequences of spin current conservation in systems with SU(2) spin symme- 
try that is spontaneously broken by partial magnetic order, using a momentum-space approach. 
The long-distance interaction is mediated by Goldstone magnons, whose interaction is expressed in 
terms of the electron Green's functions. There is also a Higgs mode, whose excitation energy can be 
calculated. The case of fast magnons obeying linear dispersion relation in three spatial dimensions 
admits nonperturbative treatment using the Gribov equation, and the solution exhibits singular 
behaviour which has an interpretation as a tower of spin-1 electronic excitations. This occurs near 
the Mott insulator state. The electrons are more free in the case of slow magnons, where the 
perturbative corrections are less singular at the thresholds. We then turn our attention to the 
problem of high- Tc superconductivity, through the discussion of the stability of the antiferromag- 
netic ground state in two spatial dimensions. We argue that this is caused by an effective mixing of 
the Goldstone (j) and Higgs h modes, which in turn is caused by an effective <j) condensation. The 
instability of the antiferromagnetic system is analyzed by studying the non-perturbative behaviour 
of the Higgs boson self-energy using the Dyson-Schwinger equations. 

PACS numbers: 11.40.-q Currents and their properties, 74.20.-z Theories and models of superconducting 
state, 74.72.-h Cuprate superconductors, 75.10.-b General theory and models of magnetic ordering 
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I. INTRODUCTION 



A. General picture and framework 

In this paper, we discuss a momentum-space, Green's function, formulation of par- 
tial magnetic ordering. Our study is motivated phenomenologically by the problem of 
high-T c superconductivity^, where the superconducting phase appears in between anti- 
ferromagnetism (AF) and the metallic phase. 

We ordinarily address magnetism (other than through phenomenological models) in terms 
of local models (e.g., the Heisenberg model) or Hamiltonians (e.g., the Hubbard model), 
which start from the picture that the electrons are almost localized at individual lattice 
sites. Their interaction, e.g., the U term in the Hubbard model, is written under this 
assumption. 

One fundamental problem which comes with this conventional approach is that when the 
magnetic order is weak and the electrons are relatively free, the assumption of locality is no 
longer viable. Specifically, the U term of the Hubbard model is written in a spin-dependent 
form in accord with the Pauli exclusion principle imposed on local wavefunctions. This is 
no longer the case when the wavefunctions are non-local. 

One would then ask, how would one formulate the origin of magnetic order, which is 
something that arises at the spatial scale of neighbouring lattice sites, in terms of something 
as delocalized as Bloch waves. The answer is that we don't. Instead of having a full 
description of the short-range effects, we parametrize them in terms of the energy difference 
AE between the spin states, and one magnetic condensate: 

<v^> M = (rM - , (i) 

i.e., the density of occupation of electrons with magnetically favourable spin, minus the 
density of occupation of electrons with magnetically unfavourable spin. (ip*ip) M is obviously 
a function of AE, and is, at least in principle, calculable for a given dispersion relation and 
total electron density. In later discussions, we shall omit the subscript g (which stands for 
'ground') and retain the subscript * (which stands for 'excited'). Thus 

(^) M = p-p^2S/a d . (2) 

S is the average spin per unit cell and a d is the volume of a unit cell. For AF, S should be 
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understood as referring to staggered spin. 

The interaction that is responsible for magnetic order is a short- distance and short-time 
effect, and is generally irrelevant when discussing long range effects. All that we require to 
know is how the condensate enters into the quantities that are relevant to the description 
of long-distance effects, in particular, the Goldstone bosons which correspond physically to 
the magnons 2 . 

The method by which we shall discuss the Goldstone boson interaction is through the 
electron Green's functions since, on one hand, p and p* are related to the integrals of the 
electron Green's functions and, on the other hand, the Green's functions are constrained^ 
by Ward-Takahashi-like identities (which, in the following, shall loosely be called the Ward- 
Takahashi identities) through the requirement of spin current conservation. 

There is also a Higgs mode, because the high-energy (high-temperature) restoration of 
spin symmetry implies the presence of a mode which becomes degenerate with the Goldstone 
bosons in the limit of zero magnetic order. 

B. The case of fast magnons in three spatial dimensions 

The case of fast (large spin- wave velocity u) magnons with linear dispersion relation (i.e., 
like anti-ferromagnetism) in three spatial dimensions is of special interest as a model case 
to illustrate our approach. 

When these conditions hold, we can write down a coupled system of Gribov equations. A 
Gribov equation^^ is a non-perturbative self-consistency equation for the electron Green's 
function, which sums the leading overlapping divergences (which plague perturbation theory) 
to all orders. 

We then solve the Gribov equations. The solution exhibits singular behaviour of the form 
(E — e)~ 3 / 2,1//2 . We interpret this as an indication of the bosonic nature of the electronic 
excitations. There is, furthermore, an infinite tower, with equal separation, of such excita- 
tions. This is a hint of the presence of a geometrical string (as opposed to the string-theory 
string which is a dynamical object with internal vibrational excitations) of variable lengths. 

We do not know how to handle the more general case of magnons with arbitrary veloc- 
ity, but the perturbative correction to the electron Green's function is not singular at the 
thresholds in the limiting case of infinitesimally small spin-wave velocity. Thus there is some 
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transition from correlated to uncorrected behaviour, depending on the value of u/v e , where 
u is the magnon velocity and v e is some typical value (~ vp) of the electron velocity. The 
limit of fast magnons corresponds to the case where the vacuum responds collectively, and 
the electronic degrees of freedom is more aptly described as being due to string-like objects, 
whereas the limit of slow magnons corresponds to the case where the vacuum does not 
respond fast enough to cause significant alteration to the behaviour of individual electrons. 

C. On high-Tc superconductivity 

Let us turn our attention to the case of two spatial dimensions, still with linear dispersion 
relation for the magnons. 

Since the infrared divergence is now stronger than in the case of three spatial dimensions, 
we expect the long-distance magnetic correlation effect to be quite radical and dominant 
when u is large. Superconductivity, in the form of electron pair condensation, thus requires 
u to be small. We demonstrate that this is usually the case, unless magnetization S is large, 
v e is small and the carrier density is small, i.e., the system is near the Mott insulator 
state. 

However, the simple exchange of magnons between electrons is unlikely to lead to pair 
formation. Specifically, magnon exchange is attractive between electrons of opposite spin. 
Thus spin singlet pairs can form, if the pair binding energy is greater than the energy that 
is required to flip the spin of one of the electrons. This spin flip goes against the magnetic 
order which is responsible for the appearance of the Goldstone magnons which mediates the 
attractive interaction, and hence pair formation is unlikely. In short, we expect that the 
pair binding energy never exceeds the spin flip energy. 

On the other hand, this restriction may be lifted for the case of an electron which is 
borrowed from a neighbouring site. As an illustration, the excitation of an e* electron with 
an up-type spin, let us say, ostensibly costs energy, but one may borrow an e electron with 
an up-type spin from a site where up- type spin is favourable, i.e., a site belonging to a 
different sub-lattice, at a lower cost in energy. 

Upon some reflection, we notice that such an effect is an instance of Goldstone-boson 
condensation. This then causes an effective Higgs-Goldstone interaction which, in turn, 
makes the AF ground state unstable against the formation of a superconducting condensate. 
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In view of this, we propose that the analysis should proceed through the computation of the 
Higgs boson Green's function to all orders. The instability should then appear in the form 
of an imaginary part of the Higgs boson Green's function. 

The behaviour of the Higgs boson Green's function is then analyzed by summing the 
finite short-distance contribution to the self-energy to all orders using a Dyson-Schwinger 
formalism. This leads to an ordinary differential equation governing the behaviour of the 
self-energy as a function of the relevant Higgs-Goldstone three-point coupling. We solve 
this equation using a power expansion and Borel summation. 

Some words of caution are necessary here. 

First, an imaginary part of a Green's function is ordinarily interpreted as leading to the 
decay of the particle. This is not the case here since there are no final state particles to 
which the Higgs boson can decay. The Higgs boson ostensibly decays into a collective state 
of Goldstone-Higgs mixture, or, the vacuum itself decays. 

Second, the instability of the vacuum does not guarantee that the stable new vacuum is 
superconducting. This statement is quite true in general, but we find that the expression 
for the imaginary part of the Green's function turns out to be of the following BCS-like^ 
form: 

A sc ~ cueute- 1 ^). (3) 

V is the 4-point interaction strength. Af(fi) is a suitable density of states. Such an exponen- 
tial form is expected in general in the case of pair-wise condensation, and hence we believe, 
even though this statement may not be rigorous, that the ground state is superconducting. 

Third, it may appear strange that the formation of superconducting pairs (which have 
charge — 2e) may be described by an imaginary part of the Higgs boson Green's function, 
when the Higgs boson is an object without electric charge. The point here is that the 
Goldstone-Higgs mixing represents the movement of electrons between the sub-lattices, and 
hence what we are summing to all orders is actually the behaviour of the electronic con- 
densate. In other words, a Higgs boson is not an exciton-like electron-hole object but the 
fluctuation in the number density of electrons. 

Fourth and last, when we discuss short-distance physics, we should in general be wary 
of the effect of short-distance interaction which, as we discussed earlier, are hidden and 
neglected. However, this applies primarily to the electrons, and not to the Goldstone and 
Higgs bosons which are collective excitations. Furthermore, the short-distance contributions 
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which we are summing to all orders are not necessarily as short as the neighbouring lattice 
sites. 

The calculated mean-magnetization dependence of the energy gap is consistent with the 
observed doping dependence of the superconducting critical temperature Tq in high-Tc com- 
pounds. Our analysis suggests the co-existence of partial AF order with superconductivity 
in the under-doped region. 



D. Organization of the paper 

This paper is organized as follows. 

We present a general discussion of spin-current conservation and the bosonic modes in 
systems with partial magnetic order in sec. [TT1 

We derive and solve the Gribov equation for electrons interacting with fast Goldstone 
magnons in three spatial dimensions in sec. IHIl 

We discuss the case of two spatial dimensions and high-Tc superconductivity in sec. [TV] 

The conclusions are stated at the end. 



II. GENERAL ANALYSIS OF SPIN CURRENT CONSERVATION 

In this section, we adapt Gribov's analysis of axial current conservation^ to the context 
of spin current conservation in systems with partial magnetic order. 

Consider the fermionic (quasi-)electron field which is an SU(2) spin doublet 

-(:)■ 

Let us define the four-vector currents. For V 1 defined by the conserved electron current 

•'electro* = , (5) 

we define 4- vector spin currents Jf (// = . . . 3, % — 1, 2, 3) of the form 

Jf = tf^f"^. (6) 

For example, J% is the up-spin current minus down-spin current. V 1 can be an operator, 
e.g., = id/dx^. 



Let us now move to the momentum space. The current vertex T M is the momentum- 
space counterpart of T M in the real space. In the spin-symmetric phase, the spin currents are 
conserved because of the conservation of electrons. That is, the following Ward-Takahashi 
identity is satisfied: 



q are d+1- vectors with components (go, q)- Qo is the energy, q is either the spatial momentum 
or the wave number (we shall be sloppy later on for the sake of the brevity of notation). 
Ai j2 refer to the spin states, but these are dummy indices here in the sense that G~ l {q, A) is 
independent of A. 

The Ward-Takahashi identity is violated in the symmetry-broken phase, since there is 
now an energy difference between the different spin states. 
We can set locally, and without loss of generality, 



ip is the spin ground state and ip* is the excited state. Even in the case of AF, we can still 
denote the spin ground state and the excited state as ip and respectively. 
Let us define the energy difference AE between the two spin states by 



If G _1 and G" 1 are both linear in energy, AE is given by e*— e and is constant up to a possible 
dependence on the spatial momentum q. AE depends on the energy qo in principle. In 
particular, at the threshold, G~ x would, in general, have a singular structure corresponding 
to the emission and absorption of the Goldstone boson <fi through the process e* — > e<fi. We 
shall encounter a case in sec. II III where the electron Green's function is strongly singular 
in the threshold regions. The electrons move in a correlated fashion, and AE cannot be 
treated as a constant under such circumstances. However, AE « const, is a reasonable 
approximation when the correlation is weak and the decay width for e* — > e<p is relatively 
small. 

After the symmetry violation, the currents Jf 2 are no longer conserved, and the Ward- 
Takahashi identity is violated by 



r"(gi - g 2 ) M oc G (qi, Ax) - G-\q 2 , A 2 ). 



(7) 




(8) 



AE = G-\q)-G; 1 (q). 



(9) 



r"(gi - ? 2 ) M oc G-\ qi , A x ) - G^(q 2 , A 2 ) ± AE (Ai ^ A 2 ). 



(10) 



s 



This AE contribution is of the same form as the coupling of the Goldstone boson 0, and 
current conservation is restored by including the contribution of the Goldstone boson. This 
is the case even when AE is not constant. Specifically, spin current conservation is restored 
for the vertex modified by the inclusion of the Goldstone boson, 



The wavy lines indicate the spin current, and the dashed line indicates the Goldstone boson. 
The two-point function in the second term is nothing but an electron loop. There is nothing 
strange in this result, since the Goldstone boson arose in the first place as the longitudinal 
component of the spin current. After taking away the longitudinal component, the remaining 
part is transverse and therefore satisfies the Ward-Takahashi identity. 
The Goldstone bosons 0i and 02 correspond to the SU(2) rotation, 



and they correspond physically to the magnons. / is the Goldstone boson form factor which, 
by virtue of eqn. (fTTT) . is calculated as the strength of the current-Goldstone-boson two-point 
amplitude. Their coupling strength with fermions is then given by f _1 AE. The Higgs boson 
mass Mh, or excitation energy, is evaluated in a similar fashion 4 to the evaluation of /. 

A. Case of quadratic dispersion relation 

For ferromagnetism, the spin-wave velocity is known to have a quadratic dispersion rela- 
tion, co oc k 2 , for small k. That is, the Green's function is given by 1 = E — ck 2 , where c 
is some constant. In this case, evaluating the electron loop in the two-point function shown 
in the second term of eqn. fTTT]) . between the spin current and the Goldstone boson, should 
yield / x (l,cq) in the soft limit such that eqn. (TIT]) satisfies the Ward-Takahashi identity. 
The form factor / has the dimension of (length) ~ d l 2 . 




(11) 




(12) 
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By considering the Oth component of the two-point function in the soft limit, we obtain 

f = -2 x J J^G{k)G.(k)AE{ky (13) 

The factor 2 is from contracting the Pauli matrices, or for the two spin orientations. The 
corrections to the electron Green's function are included in this expression. As for the 
vertex corrections, we are taking the current vertex to be bare and take the value 1 in its 
Oth component, whereas the Goldstone-boson vertex is renormalized (as otherwise there will 
be double counting). The latter correction gives rise to the renormalization of AE. f is just 
a number. This expression, and the results that follow, are therefore exact. 

By eqn. ([9]), and because the integral of a Green's function is the number density, we find 

f = 2(p-p*). (14) 

p is the number per unit volume or area. 

We can define the mean magnetization S, < S < 1/2, by 

2S/a d = p-p m , (15) 

to obtain 

f 2 = AS/a d . (16) 

Let us temporarily consider a magnetic order-parameter field $, whose rotational imparl- 
ances correspond to the Goldstone bosons. It is easy to see that / = 2v , where v is the 
vacuum expectation value of $. This implies v 2 = S/a d , which seems natural. 

B. Case of linear dispersion relation 

In the case of a linear dispersion relation for the Goldstone boson as in the AF magnons, 
let us proceed analogously to the above case of quadratic dispersion relation. Now the 
Green's function of the Goldstone boson is given by D7 1 = E 2 — h 2 u 2 \i 2 . We then see that 
the two-point function shown in the second term of eqn. (TTTT) . between the spin current and 
the Goldstone boson, is now given by fq^ (with the spin- wave velocity u = 1), where is 
the 4- momentum. In particular, the two-point function vanishes in the soft limit. In the 
case of AF, this is due to the cancellation between the two spin orientations, or between the 
two sub-lattices. 
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The form factor / 1 has the dimension of (energy xlength d ) 1//2 . Let us make use of 




-fq v wQw 



(17) 



Note that we have been rather loose in the notation here: this equation is valid only in 
the sense that the cancellation between the sub-lattices has been taken into account on the 
left-hand side. If not, we would be equating something that is finite in the soft limit on the 
left-hand side against something which vanishes on the right-hand side. 
It follows that, for small momenta, 



f 2 9^ = ~ wQw. 



For the 0-0 component of this quantity, and for nearly constant AE, we obtain 



(18) 



,9 f d d k 1 

The factor 2 is again from contracting Pauli matrices. The region of integration covers all 
points where one spin state is occupied and the other spin state is vacant. 
Let us define the mean magnetization S, < S < 1/2, again by 

2S/a d = p-p„ (20) 

We then have 

■f 2 -^ < 21 > 

In terms of the order-parameter description, / = 2v as before, and this expression seems 
natural from this perspective. 



C. Relationships among AE, (ip*ip) M , S and u 

(i/j*i/j) m is controlled by AE when AE is nearly constant. If we denote the density of 
states (of both spin states) by g(E) and if this density of states can be taken to be roughly 
constant near the chemical potential /x, we see immediately by comparing the occupation of 
the two spin states that 

p-p* = g(jj>)AE/2, (22) 
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or 

AS 

AE = "Tr^ ( 23 ) 
g(fi)a d 

Hence 

p J ^(aOA-E (quadratic dispersion relation), 
I g(/x) (linear dispersion relation). 

In sees. Ill Al and III B| we calculated f 2 through the consideration of the time component 
of the spin-current-magnon two-point function. Looking at the space component of the 
same amplitude should give us the spin-wave velocity, in principle. However, the calculation 
is messy, and it can depend on the threshold behaviour of G, G* and AE. 

A simpler and more intuitive way to obtain the spin-wave velocity, for a given form of 
the dispersion relation, is to equate the excitation energy huj^ of the ir/a magnon mode with 
the excitation energy of e to e* at n/a, because these are equivalent. This tells us, quite 
generally, that 

hu^/a, 0, 0) = e*(vr/a, 0, 0) - e^/a, 0, 0) = AE. (25) 

In particular, for a linear dispersion relation with velocity u, we obtain 

huir . ,_, aAE 4S 

= AE u = —— = - — — — -. (26) 

a nir ri7Tg(jj)a a 1 

Note that we cannot expect the dispersion relation to be exactly linear up to ir/a. This is 

an approximation. 



D. The Higgs boson 

As we expect spin symmetry to be restored in the high energy (temperature) limit, there is 
necessarily a Higgs boson h to complement the Goldstone bosons and provide the remaining 
local rotational degree of freedom. By this symmetry, the Higgs boson couples to fermions 
with the same coupling strength (f~ 1 AE) as the Goldstone bosons, but there is now an 
energy gap. This energy gap is due to corrections below the high energy scale. Explicitly, 
the corrections are given by the self-energy integral: 

/d d+1 k 
^-^(f-^AEf [G(k)G(k -q) + G*(k)G*(k - q)) . (27) 

We then compare this with the corresponding integral for the Goldstone bosons: 

^ = " 2 / (2^p^ (rlAE)2G(fc)G * (A; ~ q) - (28) 
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We see that 

E ft (g+(A£,0)) = £^( g ), (29) 

i.e., the correction to the Higgs boson at energy AE + E is the same as that to the Goldstone 
boson at energy E. Since the Goldstone bosons have zero excitation energy, the excitation 
energy for the Higgs boson is thus equal to AE (at least for almost constant AE). It follows 
that the Green's functions for the Higgs boson are given by 

, . AE (quadratic dispersion relation), 

D- h \q) = D^\ q )-{ 14 /' (30) 

I (AE) (linear dispersion relation). 

The above derivation is somewhat intuitive. Let us derive the same results more rigor- 
ously, using a UV counterterm argument. The excitation energy is ordinarily found by the 
condition that the radiative correction vanishes, i.e., £ = 0, when the excitation is on the 
mass shell. This S should be the sum of and £[/y, where £[/y is the high-energy contri- 
bution due to (if>*if)) M , and is the same between and h. The condition that + Ej/y = 
at q = completely fixes T> uv as 

ld+li 



+2 / W^H {r WGWaW = -2f- 2 AE (ri>) M . (31) 



From the requirement £/, + E^y = 0, we then obtain the excitation energy AE. That there 
is a solution to E = is a sufficient proof for the existence of the Higgs boson. 

III. THE NON-RELATIVISTIC GRIBOV EQUATION 



/ <p,q — k \ 
J » 



4>*,q ip,k ip*,q 
FIG. 1. The electron self-energy diagram. 

The exchange of Goldstone bosons is the leading long-distance contribution to the electron 
self-energy. Corresponding to fig. (TJ the self-energy is given by 

/rjd+lh. 
^^f-*AE?G(k)D{q-k), (32) 
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and similarly for S(g). The factor 2 sums over the two Goldstone modes. D is the propagator 
of the Goldstone bosons. This expression contains all corrections other than the vertex 
correction which leads to a correction to AE. 

In general, G(k) and G*(k) become singular in the threshold regions where the decays 
e* — > e<p and e — > e*<p (the latter being to a virtual e*) open up. This is fine at the 
leading order, but since corrections to G(k) and G*(k) affect AE, the general expression 
becomes troublesome, and we do not know how to handle it. This problem of 'overlapping 
divergences' is, of course, a well-known generic problem of perturbation theory. 

However, we have one exception in the case of three spatial dimensions, with a linear 
dispersion relation of the Goldstone boson, and with large u. We may, in this case, employ 
the Gribov equation framework-^. Let us therefore specialize to this situation here. 



A. Derivation of the Gribov equation 

To obtain the Gribov equation, we apply 

- - -— , (33) 
dql u 2 <9q 2 ' 

to eqn. fl32|) . The integral sign then disappears because of the following identity: 



d 2 1 d 2 



dql u 2 <9q 2 



1 lT V>( g _A0. (34) 



(Qo ~ ko) 2 — u 2 (q — k) 2 + ie u 3 

In other words, D(q) is a Green's function of the 4-dimensional Laplace equation in the 
energy-momentum space. We have assumed that the boundaries in k do not affect the 
argument, because u is large. That is, only small k is exchanged, typically. We then obtain 

d 2 1 d 2 \ (f^AE) 2 



dq$ u 2 dq 2 J * y ^ J 2n 2 h 3 u 3 

So far, this is just an algebraic manipulation which is valid at the one-loop order, but 
we now notice that something else has occurred to the right-hand side. To see this, let us 
consider the general iV-loop perturbative contribution 5E to the self-energy. The leading 
contribution is from regions with overlapping divergences, where there may be a number of 
large logarithms being involved, in the form, for example, 

~ log(n)log(r 2 ), (36) 
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where T\ and r 2 are some large ratios. Now we consider differentiating this expression twice. 
If we differentiate log(ri) once and log(r 2 ) once, we obtain something like l/r^, and both of 
the large logarithms have disappeared. On the other hand, if we differentiate only log(ri), 
the other logarithms remain intact, and so this is the leading contribution to the double 
derivative of <5£. 

Since these large logarithms are due to the integrals of D(q), it follows that the double 
derivative, when applied to the general iV-loop expression for self-energy, removes one or the 
other of the Goldstone boson propagators, and, the final contribution is of the form given in 
eqn. (1351) . The higher loop contributions amount to the vertex corrections, and so eqn. (13"5"|) 
is in fact valid to all orders in the sense of summing the leading overlapping divergences. 

Note that this approach is quite distinct from methods such as the leading-logarithm 
summation, which do not deal with the problem of overlapping divergences. More detailed 
discussion of the various contributions to the Gribov equation is available in the papers^ by 
Gribov. 

Equation fT35|) can be turned into a self-consistency equation for the electron Green's 
function in the limit of large u, since the left-hand side then becomes equivalent to the 
energy double derivative of — G~ 1 (q), i.e., 

- M G:Hq) = (J w G{q) - (37) 

Note that, in general, G^ 1 = Gq 1 — S, and we are assuming that Gq 1 has the form Gq 1 = 
go — The spatial-momentum dependence has dropped out of the expression, and we 

can write this equivalently as 

with /|j = 2n 2 h 3 u 3 f 2 . Prime refers to the energy derivative. 
B. Solution of the Gribov equation 

Let us proceed to solve eqns. ( I3"51) . We first consider the expression 



G^G; 1 )" -G^G- 1 )" = 0, 
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(39) 



which follows from eqns. (I38p . We then use the identity x 1 x" = (lnx)" + ((In a;)') 2 : 

(lnG^/G- 1 )" + (QnG; 1 )') 2 - ((hitr 1 )') 2 = 0. (40) 
Let us define z = G~ l /G~ l , so that 

(In z)" + (In z)'(ln G^G 1 )' = 0. (41) 
This can be integrated once, and we obtain 

ln(ln z)' + IniG^G- 1 ) = const., (42) 

or, 

(In z)' = GGG*. (43) 
C is a constant of integration. This implies 

G 2 = G~V, G\ = C-'z'/z 2 . (44) 

Let us make use of eqn. (1111) to eliminate the energy derivatives in eqns. (1551) . This gives 
us, almost trivially, 

G^G=(f R C)- 2 (z + z- 1 -2). (45) 

The equation for G* is obtained by substituting G* for G and z" 1 for z. 

This equation is simple, but does not seem to have a simple analytical solution. However, 
we obtain the following limiting behaviours: 

• When G diverges, z = G~ x /G~ l should also diverge. We see that the limiting expres- 
sions are given by G = (/ i? G)- 1 / 2 (-16z 3 /3) 1 / 4 and G, = (/ fl C)- 1/2 (-16/3;z) 1/4 . That 
is, G* vanishes. 

• When G* diverges, G and z vanish. The limiting expressions are given by G = 
{f R C)-^(-16z/3)^ and G, = {f R C)- l / 2 {-\Q/^fl\ 



In order to obtain a numerical solution of eqn. (1451) . it is desirable to replace z with 
something that has a limited range. In this regard, we notice that the left-hand side of 
eqn. (|45|) has a conformal symmetry with respect to transformations on z. This symmetry 
becomes more manifest when we return to eqns. (1381) and now eliminate G and G* using 
eqn. (1441) . After some elementary algebra, we obtain 

£ " i (^) 2 = + - 2 »- w 
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We then see that the left-hand side of this equation is invariant under the Mobius transfor- 
mation: 

az + b 

z — > w = (47) 

cz + d K ' 

It turns out that z is negative. Let us therefore use 

w = l^, z = -^—^- (-Kw<l). (48) 
1 - z 1 + w v--/ v/ 

We also define 



v 3 Pi = -luA-rv ( 5 °) 

dw 2 4(1 - w 2 ) 



We then obtain 



Note the similarity with eqn. ( 145]) . We may have obtained eqn. ( 15 Oj) directly from eqns. (145]) 
and (jSJ) using eqn. (148]) . 

Now y vanishes when u> = ±1, and the limiting expression is given by y — > (1 — it; 2 ) 1 / 4 . 
We should use this boundary condition to solve eqn. (150]) numerically, but this is difficult. 
We therefore use the boundary condition dy/dw = at w = and substitute some numbers 
for y at w — and, by trial and error, establish the value of y at w — that leads to the 
correct behaviour at w — ±1. The result of this calculation is shown in fig. [2] The numbers 
were obtained using the classical Runge-Kutta method with the step size of 0.0005, and 
were found to be stable against modifications to the step size. 

We found that the value of y is approximately 1.061 when w — 0. We see that y = 
(1 — w 2 ) 1 / 4 and y = 1.061(1 — w 2 ) 1 ^ are both good approximations to the actual behaviour 
of y. 



Let us define the scaled energy x as x = (E — e*)y 16/ i? 2 /3, so that y, z — > as x — > 0. 
By eqn. (T49]) . we obtain 

/Wife- (51) 

The integral on the right-hand side is finite when the integration is performed between 
w = ±1. When w = ±1, the derivative of w with respect to x vanishes, w is therefore 
periodic in x. Let us, for the sake of simplicity, adopt y = (1 — w 2 ) 1 ^. This leads to 

w = cosx, z = — tan 2 (x/2). (52) 
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y vs. \w 
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FIG. 2. y vs w. We compare the numerical results with the limiting expression y oc (1 — w 2 ) 1 ^. 
C. Nature of the solution 

We have seen in the above that the solution has the following form: 

z « -tan 2 (a;/2) (53) 

This implies that, by eqn. (1441) . 

G « i^-G- 1 tan(x/2) sec 2 (x/2), G* w =Fa/-G _1 cot(x/2)cosec 2 (x/2). (54) 
As x — )• 0, these functions behave as 

G -> i^-G- 1 ^, G* -> tV-SC -1 /^ 3 - (55) 

This is clearly a rather exotic behaviour for the Green's function of an electron. For instance, 
the integral of the Green's function, or the electron number, becomes undefinable. Note that 
this does not necessarily invalidate our starting point since we can still define (ip*ip) M as 
a finite number. Furthermore, the non-local correlator (ip*(t)ip(t + At)) is definable and is 
non-zero (at, in this case, discrete values of At). That is, the field theory is non-local. 

One interpretation of this exotic power behaviour is through the Levinson theorem which 
states that the number of states in between energies Ei and E 2 is given by the difference in 
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the phase shifts S(E 1 ,E 2 ) under scattering: 

N=-[8(E 2 )-8(E 1 )]. (56) 

71 

In this sense, we can say that there is —1/2 electron of the e type and +3/2 electron of the 
e* type at x — 0. On the other hand, when z diverges, we obtain —1/2 e* electron and +3/2 
e electron. 

This result indicates that there is no longer a physical electron. However, if we sum the 
states that are present at x — 0, we see that there is —1/2 + 3/2 = 1 state, whose charge 
is — e and spin is, if such a sum may be made, |— 1/2(1/2) + 3/2(— 1/2)| = I. In this sense, 
the electronic excitations are bosonic. 

What is the nature of these bosonic excitations? 

Our solution implies that there is a tower of such states, at each phase-space point, with 
equal separation in energy (of approximately fR7T\/3/4:) and alternating spin direction. Note 
that the number of states in each excited state is the same. This implies that the excitations 
are in the form of one-dimensional objects. This then suggests the existence of string-like 
objects (such as the trajectory drawn by spin-flipped electrons), which bind together the 
electrons. The quantized energies correspond to the length of the string. The string is not 
a dynamical object with an internal degree of freedom, as otherwise we cannot explain the 
one-dimensional nature of the excitations. 

We should add that a power-like behaviour of the Green's function is obtained also in the 
case of simple photon-exchange interaction. In this case, the Gribov equation is given by 

(G- 1 )" = ^((G- 1 )') 2 . (57) 

7T 

a ~ 1/137 is the fine structure constant. This is in the absence of screening. It is easy to 
see that this has the following general solution: 

G' 1 oc (E-e) 1/(1 ~ a/n) . (58) 

According to the above argument based on the Levinson theorem, the number of electronic 
states goes up as a result of the repulsive Coulomb interaction. There seems to be a contri- 
bution due to plasmon-like charge oscillation states here. 
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D. The case of slow magnons 



The above discussion is for the case of fast magnons with linear dispersion relation in 
three spatial dimensions. This case is special in that the leading overlapping divergences 
can be dealt with by using the Gribov equation framework. 

Let us consider relaxing one or the other of these conditions. 

First, we do not know what the situation may be with respect to changing the linear 
dispersion relation to a quadratic dispersion relation. 

Second, the infra-red divergence is stronger in 2-d. We thus expect the long-distance 
correlations to be even more dominant than in the 3-d case. On the other hand, we know 
that there is no long-distance correlation in the 1-d case because of the strong quantum 
fluctuations. Our intuitive guess (which is only a wild guess) is that the string picture is 
valid both for 3-d and 2-d, because the string is a one- dimensional object which can be 
embedded in both 3-d and 2-d. The 1-d case is obviously different in this picture. 

Third and last, let us consider the case of very slow magnons. This is interesting because 
although eqn. (I3"2"j) diverges, it is almost constant with respect to variations in the external 
4-momentum q. Explicitly, let us consider replacing eqn. (132]) by the expression 

E *Wlsmall« ~ 2 J (27r)<*+V (g - k ) 2 - u 2 k 2 + i5k - k 2 /2m ± iS { ™ ] 

in the limit of small u. k is some typical spatial momentum scale. Then this expression is 
obviously independent of q, and it can only depend mildly on go- More explicitly, we find 
that, as we would naturally expect, E* depends on g almost only through g(qo): 

^Ua^-^^l). W> 

This particular expression uses the approximation that g(E) is almost constant, e^i are the 
upper and lower limits of the electron dispersion relation, and are expected to be very far 
away from go. 

This finding, that E*, and E by a similar calculation, are almost independent of q, is quite 
intuitive. A slow magnon field tends to stay fixed at a point in space-time, and therefore 
its interaction is independent of q. When E is independent of q, there is no divergence at 
thresholds (which depend on q) and, equivalently, no singularities at thresholds in q . 

The divergent integral of eqn. ([59|) can therefore be treated as a constant to a good 
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approximation. This gives rise to the renormalization of e. AE is unaffected, and therefore 
there are no other corrections to S at any order (consider the Dyson-Schwinger equations). 

The above argument is for the limiting case of zero spin-wave velocity. This is also 
intuitive from considering the BCS theory, where the electron-pair interaction is governed 
both by the coupling constant and the cut-off at the Debye frequency ud which corresponds 
to uk in the above case. There is no superconductivity in the limit of zero ojd- What, then, 
will be the case when uk is small but non-zero? We shall argue in the next section that even 
in that case, we cannot have electron-pair formation due to magnon exchange. Therefore 
electrons remain more or less free. Note that the formation of an electron-pair condensate 
should reflect in the electron Green's function through its imaginary part (and possibly also 
the real part). 

Thus, at least for the case of linear spin-wave dispersion relation, the behaviour of the 
electrons changes qualitatively depending on whether the spin-wave velocity is fast or slow. 
For fast magnons (u 3> vp), the response of the (magnetic) vacuum is fast and collective, 
and the electrons are under its influence. For slow magnons (u <C vp), the response of the 
vacuum is slow and partial, and the electrons move around more freely. 

These results are quite intuitive, because according to eqn. u is large when S is 

large and g(p) is small. If the electron dispersion relation goes as e ~ k 2 /2m, one can see 
that the u > Vp when vf is small and S is large. In other words, the case of fast spin- wave 
velocity occurs only near the Mott insulator state or some such state with large S and small 
electron velocity. Spin- wave velocity is slow when S is small, i.e., when the magnetic order 
is broken, so that the electrons move around more freely. 

IV. A STUDY OF HIGH- Tc SUPERCONDUCTIVITY 

Since high-Tc cuprates are to a good approximation 2-d systems with AF interaction, let 
us focus on the case corresponding to these conditions. 

We saw at the end of the previous section how electrons start to become more free when 
the magnetic order is weakened. 

We then consider how superconducting pairs may be formed. 

In the exchange of a magnon between two electrons, the interaction is attractive between 
e and e#, but is repulsive between e and e, or e* and e*. This implies that for the formation 

21 



of an ee* pair, the binding energy must exceed AE, as otherwise there is insufficient energy 
to excite e to e*. But it would be unnatural that the interaction strength is so large as 
to allow such pairs to form. For instance, the BCS binding energy is expected to be less 
than the cut-off energy (the Debye frequency in BCS theory) which is given by AE. If 
such an interaction arises at all, it can only be due to something like overscreening, and 
this circumstance is more naturally associated with SDW (spin-density wave) or some such 
magnetic structure. 

We thus believe that superconductivity through magnetic exchange occurs through a 
different mechanism, and this is through, let us say, the mixing between the two (e.g., 
spin-up and spin-down) sub-lattices that constitute the AF ground state (but note that 
a completely analogous argument can be made for ferromagnetism: this argument is not 
specific to AF). 

Through this mechanism, it is possible to 'borrow' an e* electron from a neighbouring 
site to form spin-singlet pairs without paying the excitation energy AE. 

Let us estimate the size of the superconducting gap due to this mechanism, using the 
BCS expression quoted in the introduction (eqn. (J^D): 

A = 2w clrt exp(-l/W(A0), (61) 

where V is the effective 4-point interaction strength, and is estimated by f~ 2 times some 
function 0(0) of a 'mixing- angle' for the 'mixing' between the borrowed electrons and the 
initial electrons. The cut-off, cu cu t, is estimated by AE. A/"(/x) here is the number density of 
these 'borrowed' electrons. Altogether, we thus expect 

A ~ AE exp (-s(/i)/0(0)A%)) • (62) 

The expression which we shall obtain later on corresponds to Af(p) = Mh/27rh 2 , up to a 
possible constant factor. This is the number density of Higgs bosons if the chemical potential 
were exactly at the excitation energy of a k = Higgs boson. 

There is one conceptual problem with this approach, namely that this mixing between 
the sub-lattices is not an electronic excitation but essentially an alteration of the AF ground 
state, in the sense that the flipping of the spin due to this effect is associated with the 
local modification of the vacuum from, say, $ = (Q,v) to $ = (v,0), where $ is an order 
parameter. In other words, the phenomenon that we are describing here is, we believe, 

22 



Goldstone-boson condensation. What we need to discuss then is the stability of the AF 
ground state after Goldstone-boson condensation has taken place. This instability appears 
via the Higgs-Goldstone coupling which is due to the mixing of the Goldstone boson with 
the Higgs boson. This mixing is, as we shall see, a direct consequence of the Goldstone-boson 
condensation. 

A physical picture for this effect is as follows. In the AF state, moving an electron to its 
neighbouring site, without disturbing the AF order, requires flipping the spin of the electron. 
This takes place via the emission of a magnon. However, moving the electron without raising 
it to the spin-excited state corresponds to the Higgs boson, and so the Higgs boson couples 
to the Goldstone boson, or there is effectively a mixing between the Higgs boson and the 
magnon. 

We know how one might go about calculating this effect, but the full description is 
cumbersome, and we shall parametrize it by an effective mixing angle 9. This allows us to 
calculate the corrections to the Higgs boson Green's function. 

A. Higgs Goldstone coupling 

In our picture, Goldstone-boson condensation by itself costs energy, which is compen- 
sated by the superconducting condensation energy. condensation can at least formally be 
parametrized by a negative mass squared This leads to the mass matrix: 



But a Goldstone boson is necessarily massless, and so the physical modes are obtained 
by rotating and h by a mixing angle 9^ such that becomes massless. Some off-diagonal 
counter-terms are obviously necessary to stabilize the vacuum in this way, but let us skip 
the details here. 

Let us denote this angle 6^ simply as 9 in the following. 9 is determined by the require- 
ment that the total energy is minimized. That is, if superconductivity occurs, 9 is such 
that the condensation energy plus the extra energy due to exciting the Goldstone bosons is 
minimized. This is clearly too complicated, and so let us treat 9 as a free parameter in the 
following. 




(63) 
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We denote the mixed states h and 0' as 

h \ I cos 9 sin 9 . , . 

( 64 ) 

I \ — sin 9 cos 9 

In general, only one linear combination of Goldstone bosons is involved, and the other 
Goldstone boson mode remains orthogonal to both h' and </>'. 
In terms of /i^ introduced above, we have 




tantf = J-^ « ^/M h u 2 . (65) 

V M h u - v-* 

Let us assume that /i is small and the mixing angle is not huge, so that Myu 2 can be 
approximated by M^u 2 = AE as calculated in sec. Ill Dl In principle, 9 can be measured by 
measuring the discrepancy between electronic excitation energy AE and the Higgs-boson 
excitation energy, but this is probably too delicate to be useful. 
In terms of h and <j), the h<p<j) coupling is given by 

^fh(f-%<p 2 ) 2 . (66) 

We derived this from the kinetic-energy part of the effective Lagrangian density for the order 
parameter field, 



•^kinetic (^/t^i 



h 2 h 2 
2"(W = y 



u 2 (V$i 



(67) 



Of 

using f = //2 as before. Since the kinetic-energy part of the effective Lagrangian is largely 
model independent, we believe that this result is fairly general. 

We are interested in the 4> h'h coupling. From eqn. fl66|) . for near on-shell Higgs bosons, 
and in terms of h' and <fi\ the three-point coupling T is found to be given by 

^^fcXl-irfl). (68) 

The 2! is the identical particle factor. Note that this F is different from the electron-magnon 
vertex which appeared in sec. [Til Even though there is, in principle, a four-point function in 
addition to the three-point interaction, we do not know how to sum both three-point and 
four-point vertices to all orders, and therefore we neglect it here. 

One simple estimate for 9 is obtained by saying that T 2 is maximized. We find that T 2 
is maximized when 9 = n/2, which is too large to be acceptable. A local maximum occurs 
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for sin 9 = v/2/3, when 9 is 28 de grees and sin 2 #(l — 3 cos 2 9) 2 = 32/81. This seems more 
appropriate, but let us not adopt any specific value here. 

For the sake of later analysis, let us introduce the dimensionless coupling x, defined by 

-r 2 1 

X = AnMlu* X Jhuf- (69) 

Note again that this x is different from the normalized energy which appeared in sec. IIIIB1 
This is then given by 

^ J- 2 Mh>v? 2 Q . _ 2 Q)2 = AE 2 _ 2 2 ^ 

16vr(/m) 2 v ; 167r(hu) 2 g(/j,) v ' y ' 

We found in sec. Ill CI that AE = S/g(^)a 2 and that hun/a = AE. We can use these 
relations to eliminate AE and g(fi): 

sin 2 0(1 -3cos 2 fl) 2 . = -0(1(T 2 ~ lO- 1 )^. (71) 



x 



16S ' y 'S 



This is somewhat surprising in the sense that it implies that the dimensionless coupling x 
grows as the magnetization falls. However, this behaviour is correct, as we shall show later, 
since this leads to the survival of magnetic order in the large-magnetization region. 



B. One-loop results 

Our objective is to analyze vacuum instability due to the exchange of Goldstone magnons. 
In two spatial dimensions and in the leading-order approximation, the potential due to the 
single exchange of a magnon is proportional to log(r) in the real space. The potential is 
therefore unbound as r goes to infinity. This makes Higgs bosons long-distance, or infra- 
red, confined. This affects the way that the Higgs bosons can be excited. Presumably this 
disfavours single Higgs-boson excitations. 

This confinement only affects the excitations and does not affect the structure of the 
vacuum. Hence it is quite distinct from the superconductivity phenomenon which we are 
interested in, which is something that actually alters the structure of the vacuum. In short, 
we are interested in Higgs super-criticality, which occurs when the binding energy exceeds the 
energy that is needed to create a Higgs boson, and is an 0(h/MhU) = 0(a) short-distance 
effect. We should add that the super-conducting coherence length is indeed relatively short 
in the real high-T^ systems 1 , which adds support to our discussion. 
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The finite ultra-violet contributions can be summed to all orders. Let us start with the 
one-loop self-energy correction to the Higgs boson propagator. The corresponding Feynman 
diagram is shown in fig. [3j All momenta are 2 + 1- vectors, i.e., with two space and one time 
components. 



/ <p,q — k \ 

— ► ► — ► 

h, q h,k h,q 



FIG. 3. The one-loop self-energy diagram. 

In the following, let us adopt the convention h = u = 1. It is easy to reintroduce h and u 

by counting the dimensions. The square of a vector 2+1-vector is given by q 2 = q 2 — w 2 q 2 = 

q 2 — q 2 . The one- loop self-energy £i(g) is given by 

_ f d 3 k 1 1 

l{q) ~ J (2n)H (q-k) 2 + i5k 2 -M 2 + i5- (?2) 

Mq is the bare mass of the Higgs boson, and is the same quantity as that denoted as Mh> in 

the previous section. The change in the notation is so as to be more clear when discussing 

the mass renormalization. In the same way, let us call the bare coupling r . This is the 

same quantity as that denoted as F in the previous subsection. A in the above equation is 

given by 

Ao = -rg. (73) 

This is the bare coupling, to which we shall calculate the corrections which arise due to the 
vertex correction diagrams (and not due to the vacuum polarization diagrams). We have 
included the factor 2 for the two Goldstone modes. As explained below eqn. (|68|) . the factor 
2 for identical Higgs bosons has already been taken into account. The choice of the sign is 
such that positive A corresponds to repulsion and negative A to attraction. A has the 
dimension (energy) 5 x area. 

Let us say that q is small compared with the zone boundaries, so that we will not need 
to worry about the boundary effects. In this regard, u is small compared with the electron 
velocity, but it is not small here since we are not comparing it against the electron velocity. 
Since hun/a = AE, we do not have to worry about boundary effects so long as the Goldstone 
boson does not carry energy that is comparable with AE. 

26 



Eqn. (1721) can be integrated easily using standard methods^. First, the usual procedure 
of combining the denominators yields 

Z l (q) = A £ dX J ^[£ 2 + A(g 2 (l-A)-M 2 )+^]- 2 . (74) 

The Wick rotation yields 

Ex(g) = A J dX J —— [-f + X (q 2 (l - A) - M 2 ) + i6\ . (75) 



After evaluating this Euclidean integral, we obtain 

= ^ T rfA [A (M 2 - g 2 (l - A)) - iS] ^ 
Finally, under the condition M5 2 > q 2 > 0, we obtain 



(76) 



where \q\ = ^/q 2 . Near |g| = M , this diverges logarithmically, reflecting the long-distance 
confinement. In the region \q\ < M , E has the following expansion: 

S -4Sf( A «4 + 5M + "-)- < 78 > 

Let us retain the first two terms. This approximation is accurate to 10% up to \q\ /Mo = 0.7 
or so. 

Since we have the following general relation: 

G' 1 = Z-\q 2 - M 2 ) = q 2 - M 2 - E, (79) 

we can immediately write down the one-loop Z and M 2 as 

1 A M? ( A \ 

This corresponds to a (short-distance) correction to the Higgs-boson mass of the form 

Ml = Ml - E x « M 2 + = M 2 (l +x). (81) 



The dimensionless coupling x is defined in eqn. fl69|) . 

For the sake of completeness, it is relatively easy to extend the preceding analysis to 
the case with an artificial small mass m for the (pseudo-)Goldstone boson, which screens 
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the interaction with a characteristic distance K/mu. In this case, the one-loop self-energy 
correction is given by 

Ei = 7—7-7 tanh" 1 J g| (82) 
4vr|g| M + m K J 

We again see that the self-energy diverges near the threshold, corresponding to the long- 
distance correction mediated by the Goldstone boson which goes almost on shell. 

Let us now consider the vertex correction. The corresponding Feynman diagram is shown 
in fig. gj 



/ \ 
/ \ 



t| 



Mi l|r '^ h,q 2 

FIG. 4. The one-loop vertex-correction diagram. 
The one-loop amplitude is given by 

r/r--\=\[^ \ \ 1 ,o* 

1/0 °J (27r)Hk 2 + i5( qi -k) 2 -M$ + i5(q 2 -k) 2 -M$ + i5' y J 

Ti is r plus the the one-loop correction. Proceeding in a similar fashion to that for the 
self-energy, we obtain the following Feynman parameter integral: 



iyro - 1 = ^ C dx^dx^i - V a,) 

167T J Q ^ 



[(A 1 M 2 -g 1 2 A 2 A 3 ) + (A 2 M 2 -g 2 2 A 1 A 3 )-j9 2 A 1 A 2 -^] 3/2 . (84) 

In the general case, eqn. (184"|) is difficult to handle analytically, even though the first 
two integrations are relatively easy. We have been able to evaluate it with the aid of the 
on-line software^ based on Mathematica, but the result is cumbersome and so we do not 
reproduce it here. 

Here we are only interested in the constant term which arises from eqn. (184"]) . This 
corresponds to the limit in which the external momenta are small compared with M . In 
particular, when p is small, the vertex-correction amplitude is related to the self-energy 
amplitude by 

r./r.-i = ^. (85) 
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This can be deduced from the form of eqn. (1721) . We then take q 2 <C M . The application 
of eqn. (155]) to the leading term of eqn. (175]) yields 

r '/ r °- lK -ssf (86 > 

From eqn. (I80p . the one-loop estimate for the critical coupling is — A ~ 4ttMq, at which 
point the one-loop mass becomes zero, and above which the mass becomes tachyonic. At this 
value of — Ao, by eqn. (156]) . the one-loop vertex correction makes the effective vertex greater 
by 50%. Clearly, finite-order calculations are insufficient to deal with such corrections and so 
a non-perturbative method is needed. In the next subsection, we shall sum such contributions 
to all orders using the Dyson-Schwinger equations. 

Before proceeding, we note that a tachyonic Higgs mass by itself only suggests the insta- 
bility of the AF ground state, and not the form of the true ground state. For instance, one 
may very well conclude at this stage that the true ground state is given simply by v = 0, i.e., 
zero magnetization, which would be a reasonable interpretation for tachyonic Higgs mass. 
We shall show, however, that this is only partially true. In the all-order analysis, the real 
part of the Higgs-mass squared does indeed go negative above a certain value of the coupling, 
but ground-state instability occurs regardless of the sign of the real part of the Higgs-mass 
squared. 



C. Dyson Schwinger equations 

We have shown that finite-order treatment is insufficient. In particular, we would like to 
sum the finite short- distance contributions to all orders. Let us deal with this problem by 
making use of the Dyson-Schwinger equations. 

First, the self-energy diagram is shown in fig. This includes all possible perturbative 
corrections that involve the hhcf) coupling, other than vacuum polarization which we either 
discard, or assume that it is taken into account in the Goldstone boson propagator. The 
non-perturbative equation is thus a modification of eqn. ( 175]) . given by 

S = (g 2 - M 2 ) - Z- V - M 2 ) = ^ (V + £ + . . . ) . (87) 

M, E and T = y/^A are now all-order quantities. Eqn. ( 157]) sums only the finite terms in 
the higher-order corrections, and not the singular infra-red contributions. 
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h,q h,q 

FIG. 5. The self-energy diagram corresponding to the Dyson-Schwinger equation. The blobs 
represent all-order corrections. 

Let us introduce the following notation: 

m £(? = o) r 

y and w are different from the quantities that appeared in sec. IIIIBl y is the normalized 
mass, which we would like to obtain as a function of the normalized coupling x defined by 
eqn. (169|) . w is the normalized self-energy and g gives the coupling renormalization. From 
eqn. (IHTj) . we derive the following: 

w = y 2 Z- 1 - 1 = ^(y 2 - 1), xyg = w(w + l). (89) 

As for the vertex correction, we can write it as the sum of two diagrams, as shown in 

fig. EJ 



t>»p irZ rz ^p\, 



h,qi | /i,(?2 ^,91 ] \h,q 2 

FIG. 6. The two diagrams whose sum contains all of the vertex-correction contributions. 

The second of these diagrams contains a 4-point vertex, and it would ordinarily be nec- 
essary to construct another Dyson-Schwinger equation to describe it. However, as we saw 
in sec. IIV Bj. so long as we are only interested in the constant part of the corrections, the 
insertion of a Goldstone-boson vertex is formally equivalent to the derivative by Mq and 
multiplication by To- That is, the sum of the two diagrams in fig. [6] is given simply by 

30 



The partial derivative means differentiating by Mq while keeping A constant. Now S is a 
function of both A and Mq , but by counting the dimensions, we see that S must be of the 
form 

£ = M 2 w(x). (91) 

Thus eqn. f )90|) can be written as 

3x dw . . 

9 = 1 + t "-Y^' (92) 

To proceed, we make use of eqn. (1891) to eliminate (7 from eqn. (192]) . This yields the 
following ordinary, first-order and non-linear differential equation: 

3x dw . / w\ 

Y^ = ( 1 + m K 1 -^j' (93 > 

In order to obtain w = (3/4) (y 2 — 1) as a function of x, we then need to solve eqn. (193]) with 
the boundary condition 

w(0) = 0. (94) 

Unfortunately, this boundary condition makes the solution of eqn. f )93|) difficult. We do 
not think that eqn. 093 j) can be integrated analytically and, on the other hand, numerical 
implementation suffers from instability near x = where eqn. f[9"3"]) contains 0/0. 

Let us adopt a power-series solution. Simply expanding y in powers of x and substituting 
in eqn. (1931) . we are able to read off the coefficients with the boundary condition given by 
eqn. (19~4"|) . We used Maxima— for preliminary algebra and the evaluation of the first 28 
coefficients, and Form 13 after that. Using Form, we calculated the first 102 coefficients, 
which, in our implementation, is equivalent to having verified the first 101. On a Linux 
machine equipped with a 2.6 GHz Pentium 4 processor, the evaluation of the last coefficient 
took just over 5 minutes to perform. The machine has 2 GB memory, but less than 1 % of 
the memory was taken up by the calculation. The first few coefficients are given by 

y = 1 + 2 - X - l x 2 + 2 x 3 - —x A + —x 5 - 0(x 6 ). (95) 
y 3 9 81 54 K J K J 

It should be noted that x being the normalized coupling, the expansion in x should match 
with the finite terms from the perturbative calculation, order by order. For example, the 
truncation of eqn. (I95p at the first order in x reproduces the one- loop perturbative result, 
eqn. flHUl) . at the first power order in A . 
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This series is divergent for all non-zero x. In other words, the critical coupling x c , which 
is calculated as 

x c = lim - 1 , (96) 

n->oo a n 

vanishes. Here, a n stands for the coefficient of the n'th power of x. The radius of convergence 
of the series is given by \x c \. 

We plot the ratio appearing in eqn. f )96|) . for finite n, with n up to 10 in fig. [71 We have 
found that these are numerically consistent with 

-i 



a,. 



11 3n 



+ 0(n 



(97) 



. 6 2 

The coefficient of the 0(n~ l ) contribution is about 0.3, numerically. We see that for large 
n, the series is factorially divergent, and the radius of convergence is therefore zero. 
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FIG. 7. The ratio a n ^i/a n , for 1 < n < 10, compared with the expression (11/6 — 3n/2) 
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D. Borel summation 

A method which is commonly used to assign a meaning to factorially divergent expansions 
is that of Borel summation. For Y(X) given by the power series 

oo 

Y(X) = J2 A " X ~ n > ( 98 ) 

n=0 
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the Borel transformed series is given by 

oo . 

BY(t) = ^t n . (99) 

n=0 ' 

Y(X) is then recovered by the Laplace transform: 

POD 

Y(X) =A + BY(t)e- xt dt. (100) 



Since eqn. fllOOp requires X > 0, let us define 

X = l/x, A n = a n} y(x) = Y(l/x), (101) 

when x > 0, and 

X = -l/x, A n = (-l) n+1 a n , y(x) = -Y(-l/x), (102) 

when x < 0. We can then readily find the Borel transformed series BY(t). 

Alternatively, using the first two terms of eqn. ( I9~7j) . we can estimate the Borel sum to be 

BF(t)« ai (l + 3t/2)- 7 / 9 (x>0), 
BY(t) ^ ai {l-3t/2)- 7 / 9 (x<0). 

In fig. [8], we compare this estimation, with a\ — 2/3, against the series expansion using the 
actual values of a n up to aioo- For the range of BY(t) shown in fig. [SJ the contribution of 
the last term in the expansion is at most about 0.5% of the total number, and therefore 
we consider the error due to the truncation of the power series to be negligible. The two 
expressions match reasonably well. We shall henceforth make use of the estimation of 
eqns. f)103p in our discussions. 

The first of eqns. (11031) . which corresponds to positive x, has a unique Laplace transform. 
In other words, the series is Borel summable. The result of the Laplace transform is 

y(x) =a + a l (2/3) 7/9 e 2/3x x 2/9 T(2/9, 2/3x) (x > 0), (104) 

again using the approximation of taking the first two terms of eqn. (|97|) . T is the upper 
incomplete Gamma function, defined by 

POO 

r(s,a) = / t s ~ l e- l dt. (105) 

J a 
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Borel transformed series, for negative x 
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FIG. 8. BY(t) for negative x, using the actual values of a n (solid line, with n < 27), and the 
estimation of eqns. (|103p (dotted line). 

On the other hand, the second of eqns. (11031) . which corresponds to negative x, has a 
branch-point singularity at t — 2/3. In other words, the series is non-Borel summable. We 
can still derive the expression for y(x) formally and, not too surprisingly, the result turns 
out to be of the same form as eqn. (I104p . but is dependent on the choice of the contour of 
integration. Let us choose the contour such that we obtain a negative imaginary part for 
the self-energy, corresponding to the retarded Green's function. This choice of the contour 
is obtained by moving the singularity downwards in the complex t plane, as shown in fig. |9j 

IL 

•/WVVVvv* 

2/3 -iS 

FIG. 9. Moving the branch cut off the real axis. The contour of integration is along the positive 
real axis. 

This is equivalent to giving a negative and infinitesimal imaginary part to x, and so in 
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the end, we have, for both positive and negative x, as well as x = 0, 

M = l + *^(W\(w\^ r (iw h (106) 



9 \x — id J \x — id J \9' x — id J 

The definition of the incomplete Gamma function, as given by eqn. (11051) . is ambiguous for 
negative or complex a. A more general formula reads 

T(s, z) = T(s) - 7 ( S , z) = T(s) - J2 ^T^' ( 107 ) 

k=0 ' 

using which we are able to evaluate y(x) numerically. 7(5, z) is the lower incomplete Gamma 
function. The resulting form of y(x), and y 2 {x), is shown in fig. [10J 
In fig. [TUJ we also show the function 

- (4/9)r(2/9) sin(27r/9)e 2/3x, (-2/3a;)- 2/9 , (108) 

for negative x. This is an x — > —00 limit of the imaginary part of eqn. (I106p . and it agrees 
well with the imaginary part of y(x). In general, we expect that any singularity at t — to in 
the Borel transform induces a term of the form e to ^ x . There is nothing new in this divergent 
behaviour of the perturbation series, the so-called renormalon™ being a celebrated example. 
As for the exponential form of the Borel sum, since the imaginary part of the Green's function 
corresponds to the binding energy of the condensate to which the AF vacuum decays, this 
implies a BCS-like behaviour of the energy gap, as given by eqn. (RjiZl) . 

We note that, unlike in finite-order calculations, the criticality, in the sense of an imagi- 
nary part of y, develops already at x = —0, even though it is quite small until — x becomes 
comparable with 2/3. 

The real part of y remains positive even at large —x. This is against our naive expectation 
that strong attractive interaction makes the single-particle mass negative. The behaviour 
seen here suggests that there is a balance between the attractive interaction mediated by 
the Goldstone boson and the contribution due to the condensate which counteracts it. The 
former is dominant at small —x, and the latter is dominant at large — x. 

The real part of y 2 becomes negative at x ~ —0.9. Naively speaking, AF order co-exists 
with the new condensate state for x > —0.9, and is destroyed below this value of x. 
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y(x) estimated by the Borel summation technique 
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FIG. 10. y(x) (upper figure), and y 2 (x) (lower figure), evaluated as the Laplace transform of 
eqn. (fT03l) . 

E. Phenomenological discussions 

In the preceding sections, we were able to deduce the non-perturbative behaviour of the 
single-particle Green's function for the Higgs boson. 

As per the conventional wisdom^, the imaginary part of the self-energy determines the 
relaxation time of the unstable ground state which, in turn, and by the uncertainty principle, 
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is the inverse of the binding energy. The energy gap is thus given simply by 

2A = -ImM = -M u 2 lmy. (109) 

The total condensation energy, per unit area, is given by the two-dimensional phase-space 
integral 

E^L = _ f -^-ImVk 2 + M 2 = -— Im ({K 2 + M 2 ) 3//2 - M 3 ) . (110) 
A J (2n) 2 6n 

K stands for the upper limit of integration, which is close to ir/a. In the limit as the K 
tends to infinity, we obtain 

^± = -—KlmM 2 . (Ill) 

A 47r K J 

This expression, for very large —x, has the behaviour (— x) 4//9 , or proportional to A 2 . 

We do not know how to calculate Tc, since long-range order is known to disappear for 
finite temperature in two spatial dimensions. 

Let us define So by 

- 2/3x = S/S . (112) 

By eqn. (17D1. we have 



S = ^ sin 2 6(1 -3 cos 2 6) 2 } (113) 

and so S is e>(10~ 2 ) to ^(lO- 1 ). 

Using eqn. (11081) and then using eqn. (123]) . the superconducting gap is estimated as 

2A SC = 1.17AEe- s / s »(S/S y 2 / 9 = H^ e - 5 / 5o (S/S ) 7/9 . (114) 

g{fi)a 2 

This is shown in fig. [TTJ For large S, the exponential factor tends to suppress supercon- 
ductivity, whereas for small S, the AE factor suppresses it. The maximum of 2Asc as a 
function of S/So occurs at S/So = 7/9. 

Small S and large S correspond to over-doping and under-doping, respectively, and there- 
fore the general behaviour is consistent with the doping behaviour of high-Tc superconduc- 
tors, assuming that Tc is proportional to A. Out of the parameters in eqns. (11131 1114]) . 5*o 
can in principle be measured through the measurement of the magnetization at the peak 
of 2Asc, but one needs to be careful here since S refers to the magnetization when su- 
perconductivity is not present. Perhaps some extrapolation of the magnetization at higher 
temperatures can be made to the magnetization at K. This, together with the measurement 
of g(n) and a, gives us 2A SC as a function of S, to be compared against experiment. 

37 



Superconducting gap vs. staggered magnetization 
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FIG. 11. The superconducting gap, normalized to csJ 2 a 2 gF, as a function of the mean staggered 
magnetization S normalized to So- 
la the sense of the positivity of the real part of the renormalized M 2 = M^y 2 , super- 
conductivity co-exists with AF down to about x = -0.90, or S/S = 0.74 (0.74 < 7/9 = 
0.77- • •). AF disappears in the over-doping region. The pseudo-gap is naturally associated 
with the excitation of e to e„, which is characterized by AE. 

As for the overall magnitude of A sc , an 0(100 meV) superconducting gap is obtained 
when l/a 2 g(n) is (9(1 ~ 10 eV). This seems quite reasonable. The superconducting gap 
goes up when the electron density goes down (against naive expectation). That is, smaller 
values of l/m e s is favoured for the effective electron mass. 

As for the symmetry of the superconducting gap, our analysis does not single out a 
particular type of symmetry, but we expect, naively, that the requirement of forming spin- 
singlet pairs in AF background restricts pair formation in the (±7r, ±7t) directions, and so 
there are necessarily nodes in these directions. 

V. CONCLUSIONS 

We made an adaptation of Gribov's axial-current conservation techniques to the situation 
of spin-current conservation in systems with partial magnetic order. We obtained the form 
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factor, or the coupling strength of the Goldstone bosons, and the mass of the Higgs boson 
in a non-perturbative manner. We also obtained relationships among the parameters of the 
theory. Our analysis is of a general nature and rely only on conventional wisdom relating 
to quasiparticle dynamics in the momentum space, but many of the results were obtained 
using the approximation that AE, the spin excitation energy, is independent of energy or 
momentum. Some other results depend on the assumption that the inverse of the electron 
Green's function is linear in energy. These assumptions are reasonable when the spin-wave 
velocity is small. 

Our approach is distinct from methods that start from Hamiltonians with localized states 
as the bases, because the Pauli exclusion principle acting on localized states is an ingredient 
of the latter formulations. The Pauli exclusion principle is present in our approach without 
such restrictions, in the form of the Feynman propagators, i.e., in the form of the analytical 
behaviour of the Green's functions. 

We then discussed two applications of our framework. 

We first discussed the case of fast magnons, i.e., near the Mott insulating state, with 
linear dispersion relation in three spatial dimensions. This case admits a treatment using 
the Gribov equation, which sums the leading overlapping divergences to all orders. The 
solution of the Gribov equations leads to a peculiar solution in which the fermionic degrees 
of freedom seem to disappear, and we are left with something which can be interpreted as a 
tower of bosonic excitations. A natural interpretation for this phenomenon is as a string-like 
one-dimensional excitation, which may correspond to the path transgressed by a wrong-spin 
electron. This phenomenon occurs when the spin-wave velocity is large, which implies large 
magnetization and small carrier density, i.e., near the Mott insulator state. 

With increased carrier density, the system starts to lose the magnetic order, and the spin- 
wave velocity starts to fall. When the spin- wave velocity falls below the electron velocity, i.e., 
u < vf, the behaviour of the system changes qualitatively, because the collective response 
of the vacuum subsides and the electrons become more free, even though magnetization still 
remains. This cross-over occurs for reasonably large values of S. We have found that u < vf 
is satisfied unless g(n) is small, vp is small and S is large, i.e., quite near the Mott insulator 
state. This finding supports our being able to use the nearly-free electron picture to discuss 
high-Tc superconductivity. 

When discussing high-T c superconductivity, we proposed that there is a mixing between 
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the Goldstone and Higgs modes due to Goldstone-boson condensation. Goldstone-boson 
condensation is stabilized by superconductivity, and superconductivity cannot arise without 
Goldstone-boson condensation. This gives rise to an effective Higgs-Higgs-Goldstone three- 
point interaction. 

We then discussed the stability of the AF ground state with respect to the collapse due 
to this three-point interaction, which is attractive. In view of the observation of high-Tc 
superconductivity in two-dimensional AF systems, we focused on the two-dimensional case, 
which is easier in the sense that the short-distance integral converges. 

The exchange of Goldstone magnons induces non-analyticity in the Higgs-boson Green's 
function. This implies that the AF ground state is unstable at low temperature. The 
analysis of the Higgs-boson Green's function was carried out in an all-order manner, by 
using a Dyson-Schwinger summation of the relevant finite, short- distance, contribution to 
all orders. We found a solution for the mass renormalization as a power series in the 
coupling constant. The power series is factorially divergent, but can be treated using the 
Borel summation method. 

We found that an essential singularity develops in the mass renormalization when the 
coupling constant is zero. We thus conclude that the true AF ground state in two spatial 
dimensions contains a condensate into which the AF state decays. The energy gap associated 
with this condensate has the behaviour that is expected for superconductivity, and we argued 
that it indeed is the superconducting pair condensate. Superconductivity co-exists with AF 
in the under-doped region. 

A number of results of this paper can be tested experimentally. First, the expression for 
the superconducting energy gap. Second, the presence of the Higgs boson at larger values of 
carrier concentration, and the prediction of its mass. Third, the tower of bosonic states in 
AF systems at smaller values of carrier concentration, especially in the case of three spatial 
dimensions. In addition, the analytical relationships between such quantities as AE, g(/i), 
S and u are all testable. 

What might be the signals of Goldstone-boson condensation? We can at the moment 
only think of indirect signatures such as the lowering of the Higgs-boson excitation energy 
from AE, which may be too delicate to be treated as a concrete signal. 

Large T c can be achieved by making g(/i) small, or by reducing the effective mass of the 
electrons, without destroying the magnetic interaction. 
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